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ABSTRACT 



Prompt gamma-ray burst spectra evolve on short time scales, suggesting that 
time-resolved spectral fits may help diagnose the still unknown prompt emission 
mechanism. We use broad-band gamma-ray data from the RHESSI spacecraft to 
test quasi-thermal models with high signal-to-noise time-resolved spectra of nine 
bright gamma-ray bursts. In contrast to results reported in more narrow energy 
bands, the quality of the fits of quasi-thermal models is poor in relation to fits of 
the phenomenological Band function. Moreover, the best-fit parameters for the 
simplest quasi-thermal model, a black body plus a nonthermal power law, show 
significant dependence on the fit band. Models that replace the power law with 
more complicated nonthermal functions are not robust for the data considered 
here and decrease the physical relevance of the fit black body. 

Subject headings: gamma-rays: bursts 



1. Introduction 



One of the most distinct features of the initial gamma-ray emission of gamma-ray bursts 
(GRBs) is its temporal variability. Significant evolution of the lightcurve and the spectrum 
occurs on timescales shorter than the total burst duration. Accordingly, spectral fits of 
subintervals of a burst may provide improved insight into the burst emission mechanism, at 
the cost of increased complexity. 



When the iBand et al.l (1l993l ) phenomenological spectral model (the "Band function" , 
a smoothly connected broken power law) was found to successfully fit the GRB prompt 
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emission observed by BATSE, systematic time-resolved analyses of large burst samples fo- 
cused on identifying patterns in the fit para meter evolut i on in a n attempt to gain i nsight 
into the emi s sion rn echanism of long GRBs. iFord et al.l (119951 ) . ICrider et al.l (119971 ). and 
Preece et al.l (Il998bl ) considered the evolution of the peak energy Ep^ak, low-energy power- 
law index a, and high-energy power-law index /3, respectively. Broadly, these authors found 
a typical hard-to-soft decay of the emission and a general correlation of spectral hardness 
with intensity. 

The observation of hard low-energy spectral slopes {a ~ -|-1) in the initial portions of 
GRB pulses raised concerns about the viability of the synchrotron shock model, as a > —2/3 
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led a u thors to suggest such em ission could have a thermal origin (ICrider et al.lll997l : iPreece 
20001 ). iGhirlanda et al.l (120031 ) found that time- resolved BATSE LAD spectra with hard low- 
en ergy indices c ould be acceptably fit with a black-body spectrum. This result was confirmed 
by iRydd (12004J ) for hard, single pulses. In some cases the addition of a simple pow er-law 
i mpro ved the correspondence with the high-energy LAD data. These successes led iRyde 
(120051 ) to propose that all GRB emission might be decomposed into thermal and nonthermal 
components of similar magnitude, with the black-body emission providing the spectral peak 
of the prompt emission. Within the LAD band, a simple power law proved a sufficient 
approximation to the more complicated true nonthermal emission, with the resulting black- 
body plus power-law model (BB PL) fits hav i ng sim ilar values as the Band function for 
identical degrees of freedom. In iRyde et al.l (120061 ) , the authors noted that the power- law 
slope of the BBPL model avoided the OTS line of death and sought to interpret the black 
body fit re sults in the context of a photospheric model. Most recently, iRyde fc Pe'en (120081 ) 
(hereafter, |RP08| ) identified regularity in the time evolution of the temperature and fiux 
of the black-body component for single-pulse bursts and linked the no rmalization of th e 
black-body component to the size of the thermal emitting region (see also lPe'er et al.ll2007l ). 



While the simple power-law approximating the nonthermal emission in the BBPL model 
is effective over the moderate energy band (28 keV-1.8 MeV) of the BATSE LAD detec- 
tors, te sts of the BBPL mod el using data covering a broader bandpass have had mixed 
results. iMcBreen et aU kood \ successfully fit the BBPL model to GRB 041219A data from 
INTEGRAL-SPl spanning the range 20 keV-8 Me V, although t he bl ack body contributed 
only a small portion of the fiux of the main burst. iFoley et al.l (120081 ) obtained acceptable 
time-integrated fits of the BBPL model to GRBs observed by INTEGRAL-IBIS and SPI, 
although the IBIS detection require ment created a burst sa mple with relatively low peak 
spectral energies {Epeak <~ 150 keV). IGhirlanda et al.l (120071 ) considered BATSE bursts co- 



3 



observed in the X-ray by the Beppo-SAX WFC Time-resolved spectra were not available 
for the WFC data, but a summation of the time-resolved BBPL fits to the BATSE data led 
to a significant overprediction of the WFC fiux, while a similar extrapolation of the summed 
time-resolved Band function fits was much more successful. 

In this work, we have ernploye d data from the Ramaty High Energy Solar Spectroscopic 
Imager (RHESSI) (ILin et al.ll2002l ) to investigate the behavior of quasithermal spectral mod- 
els over a broad energy band. RHESSI's nine coaxial germanium detectors image the Sun at 
X-ray to gamma-ray ene rgies (3 keV-17 Me V) with excellent resolution in energy (1-5 keV) 
and time (1 binary /xs) (ISmith et al.ll2002l ). RHESSI's detectors are effectively unshielded 
above ~30 keV and receive emission from astrophysical sources like GRBs with a total effec- 
tive area of ~150 cm^. Each coaxial detector is electronically segmented into thin front and 
thick rear segments; most off-axis emission is recorded in the rear segments. The satellite 
rotates with a period of about 4 seconds. 



2. Data Analysis 

2.1. Sample Selection 

The identification of optimal fit intervals for time-resolved spectral fitting requires bal- 
ancing the need for high signal-to-noise (in order to constrain the fit parameters of the 
spectral models) with the goal o f the finest possibl e time resolution. In our previous system- 



atic analysis of RHESSI GRBs (IBellm et al.ll2008l ). we found that a S/N ratio of at least 45 



was required for adequate fitting of the most complicated models. This value i s also consis 



tent with those adopted by previous studies of time-resolved BATSE spectra ( ICrider et al. 



19971 : iKaneko et al.ll2006D . 



Because bursts may have periods of high signal which are degraded to low S/N when 
integrated over the whole burst, we chose candidate bursts from the RHESSI sampl^ with 
localizations and total S/N > 30. Radiation damage to RHESSI's germanium detectors has 
slowly degraded spectral performance; we restricted our analysis to detectors which do not 
show signs of radiation damage. We therefore considered bursts from the RHESSI launch 
in February 2002 through June 2006, after which all nine RHESSI detectors were damaged. 
Forty-five GRBs met these criteria. 

To identify subintervals within each candidate burst, we used the Bayesian Blocks algo- 
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rithm (IScargldll998l ). which identifies the most probable segmentation of a burst hghtcurve 
into intervals of constant Poisson rate. We modified the stopping criterion to generate subin- 
tervals of appropriate S/N. Since the RHESSI data are stored event by event, this "top-down" 
identification of the burst subintervals is more natural than building up data accumulated 
on fixed timescales until a minimum S/N is reached. While previous work on quasi-thermal 
models has focused on spectral evolution within pulses, our segmentation method does not 
require pulse modeling. This approach was necessary, as the bright bursts in this sample 
generally have irregular temporal structure, but it also avoids imposing additional implicit 
selection effects on the burst sample. 

We applied the Bayesian Blocks algorithm to the raw RHESSI eventlist to identify the 
most probable segmentation points. The time variation of the background, while present, 
does not dominate the GRB variability. In our modified stopping criterion, we computed 
the background-subtracted S/N (in the 60 keV-3 MeV band) for each proposed subinterval. 
If the S/N ratios of both subintervals were greater than 15, the segmentation was allowed. 
For each subinterval, we applied the Bayesian Blocks algorithm again if the S/N of that 
interval was greater than 45. After segmentation halted, we dropped any leading or trailing 
subintervals with S/N < 45 and combined any "interior" subintervals with low S/N with 
the adjacent subinterval with lower S/N. The resulting subintervals all have S/N > 45. We 
conducted spectral fitting only on bursts with at least three subintervals. 

The resulting burst sample consists of nine bursts and a total of 88 subintervals. Table 
[U lists the bursts and intervals considered in this work. Figure [T] shows the distribution 
of subinterval length and S/N. These bursts are predominantly from early in the RHESSI 
mission because for those bursts a greater number of the detectors remained unaffected by 
radiation damage. 

In order to investigate the spectral evolution of single GRB pulses in Section 14. 3[ it 
was necessary to set a lower S/N threshold. We used our modified Bayesian Blocks routine 
to identify subintervals with S/N > 25 within single, separated pulses. This "single-pulse" 
sample yielded 25 subintervals within 3 GRBs (Tabled]). 
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2.2. Spectral Fitting 



2.2.1. Fit Methods 

The RHESSI data were extracted in We fit a third-order polynomial to 

background intervals immediately before and after the burst, allowing for a potential periodic 
modulation of the background with the RHESSI rotation. We selected energy binning for 
the spectra by requiring at least ten counts in each bin of the raw source spectrum. We used 
the same energy bins for all sub intervals of each burst. Both the source and background 
spectra covered the energy range 30 keV-17 MeV. 

We determined RHESSI's spectral response by simulating monoe nergetic photons ira - 
pinging on a detailed mass model in the Monte Carlo suite MGEANT ( jSturner et al.ll2000l ). 
For a grid of angles relative to the spacecraft rotation axis, we simulated photons along 60° 
arcs in rotation angle. After interpolating to the appropriate off-axis angle, we weighted the 
60° sector responses by the burst light curve to generate the total response. Since many 
of the subintervals considered in this work were shorter than one-sixth of RHESSI's four 
second rotation period, responses for short intervals were often single sector responses. Our 
simulations indicate that the spectral response does not vary appreciably on shorter scales 
in rotation angle. 



We utilized ISIS vl.4.9 (IHouck fc Denicolall2000l ) for spectral fitting. Fitting was auto- 
mated via an ISIS script and the results stored in a database. 

We used data from the RHESSI rear segments only for all bursts except for GRBs 
021206 and 02 1008 A, which had off-axis angles small enough (18° and 50°, respectively) 
that appreciable counts were recorded in the front segments. (Coincidentally, these bursts 
also both had data decimation in the rear segments — only a fraction of the observed counts 
were recorded.) The fit energy band was 30 keV-17 MeV, with s ome exceptions: G RB 
021206 had significant atmospheric backscatter in the rear segments (jWigger et al.ll2008l ). so 
the front and rear segment data were kept separate and fit over the range 30 keV-2800 keV 
and 300 keV-17 MeV. (For all other bursts, data from all detectors were combined into a 
single spectrum.) GRBs 030519B and 040228 came through the extreme rear of the RHESSI 
cryostat (166° and 163°), where the low-energy response is more sensitive to the detailed 
mass modelling, so a fit range of 50 keV-17 MeV was employ ed. Finally, GRB 2 1008 A 
required a fit band of 100 keV-17 MeV for adequate fitting; IWigger et al.l (120081 ) found 
similar behavior in a time-integrated fit of this data. 
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2.2.2. Spectral Models 

We fit both the time-resolved and the time-integrated spectra with a variety of spectral 
models commonly used to represent the GRB prompt emission. The Band function (Band) 
i s a phenoraenolo gical model which fits the vast majority of broadband prompt spectra 



( iBand et al.lll993l ). It is a broken power law with a smooth transition between the upper 



and lower power laws: 



A{E/Epi,Y exp(-E(2 + a)/E^p,ak) E < E^^ak 

B{E/Epi^Y E > Efjreak 



with Ehreak = Eo{a - (3) and B = ^( (^z^)«-/3 exp(/3 -a). For p < -2 and a > -2, Ep^ak 
corresponds to the peak of the uE^, spectrum. The normalization A has units of photons 
cm~^ s~^ keV~^, and Epi^ is here taken to be 100 keV. 

While the Band function is often assumed to be the "true" spectral model in fits of 
GRB spectra, data obtained over a more limited energy band or at lower S/N may not be 
sufficient to constrain its parameters. We also tested other empirical models with fewer free 
parameters. The cutoff power law (CPL) is equivalent to the Band function below Ehreak'- 



Ne = A{E/E, 



piv j GXp ( 



-E{2 + a)lEp,ak) 



(2) 



We also fit a single power law (PL): A''^; = A{E / Epi^)"- . While the PL model was a poor 
fit to these spectra, its fit spectral index was useful in comparison with the other models. 

The sim plest quasi-th ermal spectral model is a black body plus power law model (BBPL). 
Proposed by iRydd ( 20051 ). it consists of Planck function and a simple power law: 



Nj. = A 



{E/E, 



piv J 



exp E/kT - 1 



B{E/Epi^y 



(3) 



The approximation of the nonthermal emission by a simple power-law i s an approxima- 
tion which is generally sufficient over th e BATSE energy band (jRvdel 120051). However, it is 
not expected to fit over a broader band ( Ryde 2004, 2005; Ryde fc Pe'ei I2OO8 ). Like the PL 
model, it diverges unphysically. 

More sophisticated quasi-thermal models replace the nonthermal PL in the BBPL with 
an empirical function with more free parameters. We attempt to fit a black body plus cutoff 
power law (BBCPL) model as well: 
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Ne = ^ ^^^E/kT-l + ^^^Z^^'™)" exp(-E(2 + a)lE,,ak) (4) 

(We retain the notation of -Epeafc for the cutoff parameter of the BBCPL model for 
clarity of comparison, but in some cases the true peak of the vFi, spectrum may be due to 
the black-body component at 3.92 kT .) However, the addition of more free fit parameters 
leads to convergence problems even with our relatively high S/N broad-band spectra. In 
addition to the unconstrained fits of the BBCPL model, then, we also perform fits in which 
the variation of the model temperature is constrained to an expected range. We describe 
the specifics of this constrained quasi-thermal model in the text below. 

Given the convergence problems faced by the BBCPL model, we have not attempted fits 
of a black body + Band model, which contains an additional free parameter. In most cases, 
the simple CPL model provides an adequate representation of the RHESSI data, with the 
Band model offering only moderate improvements. Accordingly, the BBCPL model should 
provide a sufficient test of the importance of the black-body component over the RHESSI 
band. 



3. Fit Results 

We present the time-evolution of the fit parameters of the Band and BBPL models in 
Figures [2]-[T0l The quality of the fits and the constraint of their parameters were generally 
good. For GRB 021206, the goodness of fit was poorer for some subintervals because of 
the separation of front and rear data. Slight offsets in the background subtraction created 
disagreements about the overall model normalization and increased the . We elected to 
take the averaged fluence value produced by the fit rather than introduce a fit normalization 
offset between the front and rear data for this burst. For GRB 02 1008 A, the only other 
burst with front segment data, our combination of front and rear segment spectra averaged 
out such background subtraction problems and obviated the need for such offsets. 



Wigger et al.l (120081 ) fit RHESSI data for many of these bright bursts using an identical 
mass model but independent analysis and fit procedures. Our fits of the total burst intervals 
are a good match to those results (Table [2]). In the case of GRB 021206, the presence of 
excess high-energy emission not well-modeled by the Band function increased the value 
for the fit to 17 MeV. (This component likely also contributed to the poor goodness of fit 
for this burst in the time-resolved fits, as discussed above, although its impact is harder 
to ascertain.) HETE-2 also observed three bursts in this sample. RHESSI fit higher Epeak 
values for all three bursts; in two of the cases, the peak energies were above the 400 keV 
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upper bound of the HETE energy band. 

The distributio ns of our best- fi t para meters are in general agreement with the BATSE 



results obtained by iKaneko et al.l (120061 ) . For our time-resolved Band function fits, the 
mean a is —1.00 with standard deviation 0.47 (Figure [H]). The RHESSI mean Epeak is 580 
± 280 keV, and (5 has mean —3.50 and standard deviation 0.89. The BATSE results for all 
good Band fits are a = -0.90+°j^, Ep^ak = 268^}^^ keV, and (3 = -2.33l°j^, where we have 
converted the quartile dispersions to la uncertainties under the simplifying assumption of an 
underlying normal distribution. Using Levene's test, we find that the variences of the BATSE 
and RHESSI distributions are consistent. Since the RHESSI distributions are not normal, we 
cannot compare the means using Student's t-test. The nonparametric Kolmogorov-Smirnov 
test rejects the null hypothesis BATSE and RHESSI parameter samples are drawn from the 
same distribution. However, this result is not surprising, given the different fit energy ranges 
and the few RHESSI bursts fit here. Of the three parameters, the alpha distributions are 
most consistent despite the apparent bimodality of the RHESSI sample (Fig. [1 



4. Evaluating Quasithermal Models 

We did not attempt to determine the best-fit model for each spectrum fit in this work. 
With the exception of the Band function an d the BBPL, the models fit herein have different 



numbers of fit parameters. As discussed by iProtassov et al.l (120021 ) . the conventionally used 



F-test of an additional term is inappropriate if the models being compared are not nested 
(i.e. the fit parameters of the parsimonious model are not a subset of the more complex 
model). The rigorous alternative — Bayesian hypothesis testing — would require numerical 
marginalization over our fit parameters and hence a major reworking of our analysis pro- 
cedures. Instead, we have opted to evaluate the effectiveness of the quasi-thermal models 
in three ways: first, by directly comparing goodness-of-fit for those pairs of quasi-thermal 
and empirical models with identical degrees of freedom; second, by assessing the physical 
plausibility of the model fits; and third, by considering the robustness of the fits and their 
relative consistency with each other and previous work. 



4.1. BBPL 



We begin by considering the simplest quasi-thermal model, the black body + simple 
power law ( BBPL). This model is not expected to be effective over a broad energy band (e.g., 
Rvdelboosl : iGhirlanda et al.ll2007l ). It is instructive to ask, however, how the properties of the 
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BBPL fits change when fit over a broader band. In particular, evidence of band-dependence 
calls into question the universality of results obtained even in a narrow band centered on the 
peak of the prompt emission. 

The Band and BBPL models each have four fit parameters {a, Epeak, and the nor- 
malization for Band; kT, s, and the two normalizations for the BBPL). Comparing their 
values for each subinterval (see Fig. [TTj) . we find that in 73 of 88 cases, the Band function 
has a lower value and is thus statistically favored. 

As in previous studies, we found that the peak of the blackbody emission (3.92 kT) 
in BBPL fits corresponded closely to the value of Epeak obtained in Band fits (Figure [T2b). 
Moreover, the contribution of the BBPL black-body flux to the overall flux in the band was 
significant (Figure [T5]) . 

Contrary to our expectation, we did not find systematic deviation in the residuals of the 
fit BBPL model at low or high energy. Such deviation wou ld be symptomatic of the need 



for additional spectral components (e.g.. lRyde fc Pe'eiil2008l ). 



RHESSI's broader energy band did affect the be st-fit value o f 5, the power law index of 



the BBPL model. The distribution of s reported by iRyde et al.l (120061 ) peaks around -1.5, 
while our histogram of s peaks near -1.9 (Fig. [T^ . This shift may be due to the small 
sample of bursts fit here, although as discussed in Section [3] the fit Band index a shows only 
a small shift compared to the BATSE sample. 

However, a softer fit index is expected when fitting a power law to GRB data extending 
to higher energies, as GRB spectra generically fall off more sharply above Ep^ak- The fit 
single power law index is thus a band-dependent average of the more complicated spectral 
shape. Indeed, the distribution of s obtained is quite similar to that of the index of a simple 
power-law fit to these data (Figure [T^ . and s is highly correlated with the fit PL index for 
individual intervals. Since the asymptotic spectral index must be greater than -2 below the 
peak energy and less than -2 above it, fitting a power law to a band of sufficient width which 
includes the peak energy will average to an intermediate value. 

One of the strengths of the BBPL model is that its nonthermal component typically is 
softer than the "line of death" of —2/3 predicted for optically thin synchrotron emission. 
However, this argument appears to be weakened by the sensitivity of the fit power-law index 
to the range of the data above the peak energy. It is not clear that the fit value of s provides 
useful insight into the nonthermal emission physics, as its value may be an artifact of the fit 
band. 
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4.2. BBCPL and Variants 

Given tlie cliallenges faced by tlie BBPL model in tlie RHESSI band, we tried more 
complicated nonthermal components in an attempt to find a functional broad-band qua- 
sitliermal model. We souglit to evaluate tfie efficacy and robustness of such models and to 
assess whether the thermal component retained physical relevance and realism. 

After a simple power-law, the cutoff-power law (GPL) model is the next most basic 
nonthermal model. Unlike the PL model, it has the advantage of providing a generally 
adequate representation of the RHESSI data considered here over the full band. Accordingly, 
we tested a black-body plus cutoff power law (BBGPL) model. 

Given the success of the basic GPL model at fitting these data, it is perhaps not surpris- 
ing that the GPL component of the BBGPL model dominates the fits. Figure [T2b shows that 
the Epeak of the BBGPL fits closely corresponds to the peak energy fit by the Band func- 
tion. The peak energy of the black-body component alone, however, is almost completely 
uncorrelated with the peak energy of the phenomenological models (Figures [T2li-f). This 
behavior is in marked contrast to the success of the BBPL model of identifying the spectral 
peak with the black-body peak. Moreover, many of the fit values of kT in the BBGPL model 
are limits only, indicating the lack of robustness in the fits. 

The bolometric black-body flux is obtained by integrating the differential energy spec- 
trum of the Planck function ENe over all energies; it equals A{kT)^Ep^^'ir^/15. Its value 
is fairly similar for the BBPL and BBGPL models. However, since the best-fit black-body 
temperature is often near the limits of the RHESSI band, the fraction of the flux in the 
RHESSI band provided by the black body drops sharply in the BBGPL model (Figure [T5]l . 
The black-body component is accordingly of less relevance in explaining the observed data. 



Following iGhirlanda et al.l (120071 ). we also attempted "constrained" fits of the BBGPL 
model in an attempt to force the black-body component to provide the peak of the spectral 
emission. We fit a BBGPL model with the black-body peak energy (3.92 kT) constrained 
to lie within the 90% confidence limit of the best-fit Band Epeak valued As expected, the 
value for the constrained BBGPL fits increases dramatically relative to the BBGPL fits, 
and the best-fit kT value typically does not have 90% confidence limits within the allowed 
band. The black-body contribution to the total flux in the RHESSI band is not significantly 
different than for the regular BBGPL model and remains small. 



■^Our approach is less restrictive than that of 



Ghirlanda et al. 



( 20071 ). who fixed the black-body peak at 
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4.3. Time Evolution 



One of the most notable features of the BATSE BBPL fits presented by IRP08I is the 
remarkably universal behavior of the time evolution obtained in the BBPL fits. For their 
sample of single-pulse bursts, they found that both the fit black-body temperature and the 
bolometric black-body flux exhibit a regular power law rise and decay. They found identical 
break times for the flux and temperature evolution in a given burst. Moreover, a dimension- 
less proxy for the black-body normalization, 7^(^) = {Fbb / (rT'^f'^ = 1.01 x IQ-^^A^/^E'^^, 
is related to the size of the photospheric emission region and increases monotonically even 
for bursts with complex light curves. 

We consider the time evolution of the black-body components o f the R HESSI BBPL 



and BBCPL fits in an attempt to reproduce these behaviors. Since the IRP08I time-evolution 
results were obtained from bursts with simple single-pulse time profiles, we used a single-pulse 
burst sample with slightly lower signal-to-noise ratio than our high S/N sample of mostly 
complex bursts. Unfortunately, only 3 bursts met the selection criteria described in Section 
O These bursts, GRBs 020715, 030329A, and 060805B, have 8, 12, and 5 subintervals, 
respectively. Accordingly, caveats of small sample size and limited time resolution apply to 
our analysis of their time evolution. 

We plot the time evolution of the temperature, bolometric black-body flux, and TZit) 
in log-log space in Figures [15], [16], and [T7J Within the moderate time resolution of these 
data, there is a clear suggestion of the correlated power-law rise, break, and decay of the 



temperature and black-body flux for the BBPL model, as found by IRPOS 



For these three bursts, th e BBP L TZ{t) values are generally consistent with a constant 
value. While a minority of the RPOsI sample showed such behavior, far more frequently the 



RPOSi bursts had 7^(t) values which increased as a power-law over an order of magnitude. 



The authors interp reted this increase in terms of the expansion of the photospheric emis- 



sion region. 



RP08I also found frequent monotonic increases in TZ for bursts with complex 
lightcurves. Our high S/N sample, consisting primarily of bursts with such lightcurves, also 
all have TZ (t) consistent with a constant. The values of TZ we obtain are consistent with 



the RPOsI sample, however. Both the single-peak and the high S/N samples range from 



7 X 10-2°-2 X 10-18 and peak near 3 x 10"^^ 

The time evolution of the BBCPL black body exhibits few clear trends, as expected 
given the black body's convergence problems and minimal contribution to the fit. For the 
single-pulse sample, the bolometric black-body flux seems to track that obtained in the 
BBPL model (Figures [T514T71) . There is less correlation between the temperatures of the two 
models, though, and the time evolution of TZit) for both the single-pulse sample and the 
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high S/N sample shows httle temporal regularity. Moreover, the BBCPL TZ(t) value can 
be more than an order of magnitude different than that obtained in the BBPL fits, ranging 
from 10~^^-10~^^. Such changes would imply orders of magnitude change in the inferred 
size of photospheric emission region, casting further doubt on the physical significance of the 
black-body component of our BBCPL fits. 



5. Conclusions 

We have tested quasi-thermal spectral models of GRB prompt emission over a broader 
energy band (30 keV-17 MeV) than previously considered. Many of the successes of the 
BBPL model in the BATSE band (acceptable xl, physical spectral slopes for the nonthermal 
component, universal time evolution) are not reproduced consistently in the RHESSI data 
considered here. We have attempted to construct more realistic quasi-thermal models by 
combining a black body with a cutoff power-law (BBCPL). However, even using our intervals 
selected for high quality (S/N> 45), the BBCPL model exhibits convergence difficulties. 
More crucially, though, the fit black-body component loses its relevance as the origin of the 
spectral peak. It is possible that bursts with simpler light curves might exhibit better fit 
behavior or that data obtained over an even wider energy band could clarify the transitions 
between thermal and nonthermal emission. Data from Fermi-GBM and LAT will provide 
useful tests. However, results from GRB 080916C, the first bright long burst observed by 
both instruments, indicated that a Band function pro vided good repres entation of the burst 



spectrum over the full 8 keV-200 GeV energy range (lAbdo et al.l 120091 ). 



The challenges of untangling the origins of the prompt GRB emission from the prompt 
gamma-rays alone are long-standing. New observations of higher energy gamma-rays by 
Fermi-LAT , of simultaneous long-wavelength emission by transient monitoring campaigns, 
and of gamma-ray polarization by Compton telescopes will enrich our view of the prompt 
emission and provide new clues to its origin. 
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Table 1. Intervals used for spectral fitting. The intervals are quoted relative to the start time Tq. 



GRB Off- Axis Angle Tq (UT) Interval Edges (sec) 



Ili-li S/X Samplo 


020715 


72° 


19:21:08.666 


0.000, 0.356, 0.646, 3.336 


021008A 


50° 


07:01:01.399 


0.000, 1.265, 1.985, 2.541, 2.727, 2.934, 3.142, 3.278, 








3.542, 3.889, 4.054, 4.314, 4.438, 4.995, 8.153, 8.617 


021206 


18° 


22:49:12.951 


0.000, 1.983, 2.163, 2.267, 2.405, 2.626, 3.012, 3.176, 








3.377, 3.487, 3.559, 3.629, 3.714, 3.809, 3.871, 3.959, 








4.024, 4.099, 4.361, 4.618, 4.699, 4.792, 4.945, 5.068, 








5.168, 5.346, 5.552, 5.718, 6.112, 6.217, 6.603, 6.906 


030329A 


144° 


11:37:25.954 


0.000, 2.266, 2.833, 3.520, 3.979, 4.351, 5.019, 6.226, 








13.431, 14.335, 15.006, 15.687, 17.553 


030519B 


166° 


14:04:52.951 


0.000, 5.769, 6.352, 6.726, 7.027, 7.214, 7.625, 11.065 


031027 


102° 


17:07:37.744 


0.000, 3.236, 4.774, 6.065, 9.469 


031111 


156° 


16:45:20.776 


0.000, 0.329, 0.971, 1.329 


040228 


163° 


00:08:50.910 


0.000, 4.751, 6.291, 8.458, 13.995, 15.009, 17.221, 








17.563, 18.338, 22.981, 27.748 


040810 


144° 


14:16:40.241 


0.000, 17.055, 24.995, 45.637 


Single-Pulse Sample 


020715 


72° 


19:21:07.907 


0.000, 0.543, 0.759, 1.115, 1.236, 1.397, 1.521, 1.810, 








5.135 


030329A 


144° 


11:37:25.954 


0.000, 2.280, 2.895, 3.356, 3.997, 4.355, 5.020, 5.287, 








5.743, 6.287, 7.529, 9.616, 13.523 


060805B 


97° 


14:27:15.989 


0.000, 1.245, 1.632, 2.537, 3.428, 4.451 
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Table 2. Comparison of time-integrated Band fit parameters for bursts fit in other works. Results of the closest 
analagous fit are shown. Only the limits of the RHESSI fit energy band are presented when different ranges were used 
for the front and rear segments. Parameter uncertainties have been converted to 90% C.L. where necessary. 



GRB Instrument Fit Band a. EpeakO^eV) P X^/dof Reference 



020715 


RHESSI 


30 kcV-17 McV 


— U. iU 


477+1 


-2 79+0-24 

— U.oU 


34.5/26 


(1) 




RHESSI 


30 keV-15.7 MeV 


-0.776±0.072 


531±39 


-3.14±0.41 


106.3/113 


(2) 


021008A 


RHESSI 


100 keV-17 MeV 


-1.360+n34 


733+i^ 

— ID 


~3 . 59 A 1 /I 


59.7/30 


(1) 




RHESSI 


300 keV-15.7 MeV 


-1.493±0.092 


677±54 


-3.73±0.30 


79.5/84 


(2) 


021206 


RHESSI 


30 keV-17 MeV 


-0.817+niJ 


718+i? 


-2.912+ni 


118/68 


(1) 




RHESSI 


70 keV-4.5 MeV 


-0.692±0.033 


711±12 


-3.19±0.07 


176.5/174 


(2) 


030329A 


RHESSI 


30 keV-17 MeV 


-1 648+°-°^° 


1391^ 


-3 03+°-2^ 


101.4/46 


(1) 




RHESSI (peak 1) 

V J- / 


34 keV-10 MeV 


-1.608±0.063 


157.2±8.6 


-3.48±0.87 


84.3/90 


(2) 




RHESSI (peak 2) 


34 keV-7 MeV 


-1.78±0.11 


85±18 


-3.04±0.49 


103.3/83 

1 


(2) 




HETE 


2 keV-400 keV 


-1.32±0.02 


70.2±2.3 


-2.44±0.08 




(3) 




HETE 


2 keV-400 keV 


-1 26+°°^ 

-'-•^'^-0.02 


68±2 


9 00+0.05 
-Z.ZO_o.o6 


213.7/139 


(4) 


030519B 


RHESSI 


50 keV-17 MeV 


1 171+0.049 
-1.1 / l_o.o46 


43711 


q 1 1 +0.23 
-o.ii_o,4o 


66/35 


(1) 




RHESSI 


70 keV-15 MeV 


-1.048±0.069 


417.2±21 


-3.11±0.30 


86.3/75 


(2) 




HETE 


2 koV 100 keV 


-0.8±0T 




-1.7±0.2 


92.0/121 


(1) 


031027 


RHESSI 


30 keV-17 MeV 


-u./oy_o.o80 


31911 


< -2.79 


63.9/56 


(1) 




RHESSI 


60 keV-6 MeV 


-0.94±0.13 


338±25 






(2) 


031111 


RHESSI 


30 keV-17 MeV 


^■'^ ' -"--0.061 


8411^^ 


o Qq+0.22 


50.2/56 


(1) 




RHESSI 


38 keV-15.7 MeV 


-1.102±0.059 


844±97 


-2.364±0.11 


128.3/117 


(2) 




HETE 


2 keV-400 keV 


-0.82l°:f 


404+^? 






(5) 


040228 


RHESSI 


50 keV-17 MeV 


1 «n+0.037 
-'-•"'J-0.034 


7691^^ 


^•^•^-0.36 


107.9/84 


(1) 


040810 


RHESSI 


30 keV-17 MeV 


1 4c:+0.12 


3211^^ 


< -2.52 


55.7/56 


(1) 



Note. 



_m This work (2) IWigger et al.l (|2008|) (3) IVanderspek et all (|200J) (4) ISakamoto et al.l (120051 ) 



Pelangeon et al.l (120081 ) 
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Fig. 1. — Distribution of subinterval length and signal-to-noise ratio for the high S/N burst 
sample considered in this work. 
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Fig. 2.— Spectral evolution of GRB 020715. The 60 keV-3 MeV lightcurve has 0.05 s time 
bins; counts are shown only for those segments used in the fits. The fit model parameters 
are black for the Band model and gray for the BBPL model. All errors are 90% C.L. Peak 
energies are Ep^ak for the Band model and 3.92 kT for the BBPL model. The Band model a 
is plotted with the power-law index s of the BBPL model. The errors on the fraction of the 
flux provided by the black body in the RHESSI band were estimated from a Monte Carlo 
bootstrap. Points where the value was unconstrained in the fit are omitted from the plot. 
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Fig. 3.— Spectral evolution of GRB 021008A. Symbols as in Fig. d 
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Spectral evolution of GRB 021206. Symbols as in Fig. [2l 
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Fig. 5.— 



Spectral evolution of GRB 030329. Symbols as in Fig. [2l 
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Fig. 6.— 



Spectral evolution of GRB 030519. Symbols as in Fig. [2l 
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Fig. 7. — Spectral evolution of GRB 031027. Symbols as in Fig. [21 The lightcurve appears 
flat because the background level is not visible in the time interval plotted. 
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Fig. 8. — Spectral evolution of GRB 031111. Symbols as in Fig. [2l 



- 28 - 



1.5x10 




i_ 

C 
LJJ 

CO 




X 
CD 
T3 
C 

i 



o 



■a 
c 

CO 

m 



1500 r 

1000 7 

500 E" 


-1.0 
-1.2 

-1.4 

-1.6 

-1.8 

-2.0 

-2.2 
-2.4 
-1.5 

-2.0 
-2.5 
-3.0 F 
-3.5 '- 

-4.0 §- 

-4.5 
1.0 



o 



m 
m 



0.8 
0.6 
0.4 

0.2 b 

0.0 
1.6 

1.4 b 

1.2 

1.0 

0.8 
0.6 



--+ 

1— 



1+ 



-+_ 



Jl- 



tt" 



tt 



□ 



5 10 15 20 

Time (sec) 



25 30 



Fig. 9.— Spectral evolution of GRB 040228. Symbols as in Fig. H 
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Fig. 10.— Spectral evolution of GRB 040810. Symbols as in Fig. H 
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Fig. 11. — Comparison of the reduced chi-squared values for fits to the Band and BBPL 
models, which have identical degrees of freedom. The Band model is generally preferred. 
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Band Epeak (keV) BBCPL kT * 3.92 (keV) 

Fig. 12. — Comparison of peak energies obtained in Band, CPL, BBPL, and BBCPL fits. 
As expected, there is good agreement between the Epeak values of the Band and CPL models 
(panel a). Notably, the black-body peak (3.92 kT) of the BBPL model and the BBCPL 
Epeak parameter also match the Band function Epeak (panels b and c). However, the black- 
body peak of the BBCPL model shows little correlation with the spectral peak identified by 
the other models (panels d-f). A few points which were completely unconstrained in one 
dimension are omitted in the plot. 
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BB Flux Fraction 

Fig. 13. — Histogram of the proportion of flux in the RHESSI band contributed by the black- 
body component in the BBPL (black) and BBCPL (gray) models. With the nonthermal 
component providing the spectral break in the BBCPL model, the relevance of the black 
body component diminishes. 
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Fig. 14. — Histogram of the fit low-energy power law indices for the Band model (solid 
black), BBPL model (solid gray), and simple power- law model (dashed black). The indices 
of the PL model and the BBPL model are highly correlated. The apparent bimodality of 
the Band a index is likely an artifact of the limited burst sample (Section [3]). 
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Fig. 15. — Time evolution of the temperature, black-body flux, and TZ for the single- pulse 
burst GRB 020715. BBPL model fits are black, while BBCPL fits are gray. 
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Fig. 16. — Data for the single-pulse burst GRB 030329. Symbols as in Fig. [151 
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Fig. 17. — Data for the single-pulse burst GRB 060805. Symbols as in Fig. [151 



